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SUMMARY 


A procedure is presented for obtaining "idle pressure distri- 
bution in a two-dimensional, incompressible, and nonviscous flcfw 
on an arbitrary airfoil section In cascade. The method considers 
directly the influence on a given airfoil of the rest of the 
cascade and evaluates this interference by an iterative process, 
which appeared to converge rapidly in the cases tried (about unit 
solidity, stagger angles of 0 ° and . Two variations of the 
basic interference calculations are described. One, which is 
accurate enoti^ for most purposes, involves the substitution of 
sources, sinks, and vortices for the interfering airfoils; the 
other, which may be desirable for the final approximation, 
involves a contour integration. The computations are simplified 
by the use of a chart presented by Betz in a related paper. EEhe 
numerical labor involved, while considerable, is less than that 
required by the present methods of conformal transformation. 
Illustrative examples are included. 


INTRODUCTION 


The rapid Increase of interest in the design of fans and 
turbines has led to many studies of the tTO-dimensional flow past 
infinite lattices. Most of these studies involve approximate 
procedures (for example, references 1 to 3) or present solutions 
for special classes of shapes (references it- and 5) • Recently, 
attempts have been made to obtain exact solutions by conformal 
transformation of the lattice to a circle. To this end, Howell 
(reference 6 ) used a procedure that first transformed the lattice 
to an isolated S-shape figure, which could then be transformed to 
a near circle by successive Joukowski transformations and finally 
to a circle by the method of reference 7* In reference 8 the 
cascade was transformed first to a near circle an6. then to a circle, 
also with the use of several stages of conformal mapping. In 
reference 9 the lattice was mapped into a lattice of straight 
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parallel , lines by means of a function -tiaat was determined, with iSie 
aid of the transformation of this line lattice to a circle. (See 
references 10 and 11.) These transformations are of considerable 
interest, theo3:*eticaliy . The methods of references 6 and 8 require 
lengthy coa^utations , however, and difficulty has been experienced 
in obtaining accurate mmierical results with the- method of 
reference 9. All three methods require modifications for hi^ly 
cambered contours or for lattices of high stagger and solidity. 

The method presented heredJi does not seek a conformal 
transformation directly but, like ihe older approximate methods, 
seeks to evaluate the interference at each airfoil due to the 
presence of all the other airfoils of the cascade. The velocity 
distribution on each airfoil is considered to be the sum of that 
corresponding to its presence in the uniform free-stream flow plua 
that corresponding to its presence in the interference flow. The 
tnterfereiice 1s calculated' from ihe velocity distribution on the 
airfoils so that the method reduces to on itei’ation process in ■rfilch, 
for the first approximation, the interference is computed by assuming 
the .free-str.eam velocity distribution to exist on each airfoil, and 
in subsequent approximations' this velocity is corrected according to 
the interference derived in the preceding approximation. A solution 
is thus found for an arbitrarily specified angle of attack, and this 
solution is used to find the conformal transformation to the circle 
and thence the solution for any other angle of attack. 

The present .method has been fotand appreciably less laborious 
than the methods that seek the conformal transformation directly 
and is also considered more flexible in that it may be adapted to 
a variety of cascade problems that would be difficult to solve by 
formal transformation methods; for example, the problem of the flow 
about double cascades (or superimposed lattices) or certain types 
of inverse problems involving the determination of the sotting or 
solidity for a given airfoil in cascade. Some of the features of tlie 
Interference .and Iteration methods used should also be useful in the 
solution of flows involving a finite number of interfering bodies. 


S'XMBOIiB 


V flow function (complex potential) 

velocity potential 

$ 


stream function 
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V velocity at izifinlty 

r circulation 

K mapping-function parameter 

V local velocity 

7 vortex strength 

m source strength 

z complex variable of physical plane (x + iy) 

z* fixed point In physical plans 

$ complex variable of reference plane (^ + iT) 

t 

c profile chord 

c^ profile chord used ip transformation of refei'snce 7 

d cascade spacing (distance between corresponding points on 

adjacent blades^ see fig. 1) 

9 central angle of perfect circle obtained in transformation 

of reference 7 

0 central angle of unit circle of figure 1 

s surface length on profile 

P blade angle (angle between stagger line and normal to 

chords; see fig. l) 

CT solidity (ratio of chord to distance between profiles) 

X angle between flow direction and normal to stagger line 

a angle of attack relative to blade chord 

T angle of zero lift for cascade, relative to blade chord 

Ap static pressure rise 

CO turning angle of flow 
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p density of fluid 

Fourier series coefficients 

Subscripts; 
f free stream 

d disturbance 

c compensating 

r due to circulation change 

a additional 

T total 

t tall stagnation point 

n nose stagnation point 

T.E, trailing edge 

B due to source rows 

V due to Tortex roira 

z physical plane 

^ reference plane 

0 mean flow 

1 incoming flow .. 

2 outgoing flow .. 

Xq at flow direction X^ 

Xq* at flow direction Xq* 

THEORy OP INTEEFERENCE CALCULATIONS 

In order to explain better the basic concepts and procedures of 
the interference calculations, discussion of the iteration steps will 
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■be postponed for the present, and the interference calculations will 
he described as if thej were being used to verity a known solution. 

Breakdown of the flow function into four ooinrponentB . - Attention 
is fixed on one airfoil of the infinite cascade which will be 
designated the central airfoil. The flow function on the boundary 
of this airfoil is considered to be the sum. of the following components: 

Vf the flow function for the central airfoil, considered as isolated 
in the free -stream flow (the vector average of the flow far 
in front of the cascade and the flbw far behind the cascade) • 
Inasmuch as the boundary is a streamline in this flow, 

Wf = . 

the disturbance along the contour caused by the presence of 
all the othez’ airfoils of the cascade, designated the 

external alrfoy .8 

the compensating flow function (which may have singularities 
only within the central airfoil) that is req.uired to 
maintain the airfoil a streamline in the presence of the 
dlstxirbance' flow. It is determined by the condition that, 
on the boundary, its stream fmction must be eq.ual and 
opposite to the disturbance stream function* Thizs, 

Wc = ®o + ■ vhere , 

Wp the contribution of the circulation that must be added to 

Bialntain the. tr ailing-edge condition; it has only a real 
component 

The Bim WjL + + ¥p represents the net change of flow 

function due to the presence of the external airfoils; it will be 
designated the additional flow function ¥a = ^f 

will be designated the total flow function Wq, = 

The evaluation of the isolated, or free-stream, flow is 

readily performed by the method of reference 7 and req.uires no 
further discTisslon in the present paper. The disturbance flow can 
be calcTzlated when the potential distribution (or velocity 
distribution) on the external airfoils is known. Finally, the 
compensating flow and the circulation flow are readily determined, 
as will be shown, when the disturbance flow is known. In the 
following sections two methods of calculating the disturbance flow 
will be described: the approximate source-vortex method and the 

exact contour-integral method. 





6 


MCA TN JTo. 1252 


Plsturliance flow 'by ap-proxliaata soiirce -vortex method . ~ Each 
of the externa,! airfoils is considerecL to be,.. adequately represented 
by an .arrangement of about 'two sources three slides (or negative 
aoTurces), and f Ive ■vortle'es distributed along its mean line. The 
strengths and locations of these singulari1^.es are: .chosen on the 
. basis of sh.o^^i'wis.e'.-ihickness distributlb:n\ah3.~pha^^ 
velocity distributioii.^ ‘The choice is somei&at, arbitrary and may 
be left to .the Judge^ht of the worker; however, a detailed method 
_,of choice. has been' 'd^hci*ibed in the section entitled "Ccmputatlonal 
Methods." The disturbance flow, then, is..-&at.of about ten 
infinite rows' of >. singularities, equally spaced along the cascade 
direction except -that none are located wiiere the .central airfoil is 
'to be placed. The field’bf each vortex row is shown,' in figure 2 
vbere, for corivenience, -tiie vortices are assumed -to :be,v of unit 
streng-th, spaced at unit die-bance alcaig the y-axis. This figure is 
from reference 1 and the eq^iation for the floWj i^_ (reference 2) 

- ... ‘ -V ■ - - ^ .i-.- * f 

*1 N • i'v’.'J'O • J 

\f a logg sinh Jtz + .1- logg.'Jtz 

•. ,2jt , 2it 

In order .'bo. find, the contribution to tho; distj^ba^ flow 
caused by a-, row of vortices at, say, .0.3 . 9 hor^'on .the'ex'bernal air- 
foils , -bha . .pentral airfoil j ■■■drawii 'bo scale and’ pro^Aj . or ien-bod 
relative to the 'cascade dlrectiohy- is pl^e'd At cen^r of 
figure 2, .with ‘the origin at 0.3<Jhord on'tji'e mean'.lln^ The values 
of velocity potential, and stre'Sm f^uhctlon ?e'e^d a.t seleclwd points 
eilong ’liie airfoil contour ,.'biultiplied by the nssumfed vortex strongth, 
give directly the contribution of -this vortex to. aivd 

By shifting tlie oen.'bral airfoil so. ^at iiie origin is located, in 
tinn, at. each- of the- other assumed vor-tex positions, along the moan 
line and arepeatlng -the foregoing process, the con-tributiaas of all 
'the vortices in -the external airfoils are ob'talned at the same 
points. ...The suiii of these values at a giVen point c>n the ■ cen'feral 
airfoil represents ;the 'contribution of _'th© ;yortex singi^ in 

•the lattice -bo the dis'burbance fmctioh at ’■that point . The 

contributions of..- the sources can be found in -the same way except 
•that the lines marked ' are considered as -$ and "tlie lines marked 
$ are considered, as ^ • Sinks are considered as negative sources. 

Contour -integral method for evaluatlnfi distur'bance flow _ 
function. - ‘In the preceding. section/ the -'dlsturbahce f ielcL 
calcula-bed appi^Xitia;^’;!^ .by representing •'eShh airfo'ii by 'a somewhat 
arbl'brary arrang^ent .of vortices, sources, and sinks dis'tribu'bed 
on the mean lihe,.. .An .airfoil may be represented exactly by a 
continuous .dis-ferlbution of vertices along its contour, the lineeir 
density cf which at .every point equals the velocity on the airfoil 
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at that point (reference 12). 0?he field at a point on the central 
airfoil due to a row of corresponding surface elements of the 
external airfoils (that is, a row of vortices of strength v^i ds) 

may he obtained directly frcm figure 2. Integration of this 
contribution along the contours of the external airfoils provides 
an exact deteirminatlon of the disturbance field. The procedui‘6 is 
an obvious modification of the preceding approximate method. 


Let ® and $ (without subscripts) denote, respectively, the 
potential and stream function, of the row of unit vortices in 
figure 2. In ordei- to determine the disturbance potential and 
stream, function at a point z' on the central airfoil, the air- 
foil contour, drawn to scale and correctly oriented relative to 
the cascade direction, is superimposed on figm*e 2 so that the 
origin falls, in turn, at a number of points z on the contour, 
and for each setting values of 5 and f are read at the 
point z ' . Then the disturbance flow function at z ' is given by 


<&^(z') ^ J ^ Tip(z) ds 





ds 



where 

V|j,(z) local velocity on the, airfoil at variable point z 
s distance along airfoil cbntour 

$ valTses read at z ' when origin of figure 2 is at z 


and the Integration is performed along' the, airfoil contour. 

Since Vy(z) ds = d®ip(z), the foregoing, equations can be 

rewritten as . ■ , . • - „ 

®a.(z') -\fj^ 

■ ?ji(z')- d®y(z) 

so that the disturbance potential and stream function at point z' 
•are readily evaluated by plotting. and : $ ■ -against $rji and 

measuring the area under the curves. 
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D etermination of ocanoensatlng flov and circulation flow .- As % 

has been indicated, the compensating flow function may have 
singularities only within the central airfoil contour, .end on the 
contour, the stream function must be exactly equal and opposite to 
the disturbance stream function. From the known transformation 
of the isolated airfoil to the circle, which was found in the 
process of determining W^, the correspondence between points on 
the airfoil and points on the circle is known. If, then, the 
desired. compensating stream fuocticai is expanded as a Fourier series 
in terms of the circle angle, cp. 


H 


CO 

^ (a^ cos ixp + b^ sin rxp) 

its corresponding velocity potential will be (reference 7) 
conjugate series 

c» 

Z_ (-b^ cos nq> + sin ikp) 

The determination of from is readily acccmpllshed by tha 
method of reference 13 . 

In order to maintain the trailing-edge condition, a vortex rg_ 

must be added at the center of the circle of such strength that 

rg,/2« equals the value of -d^ /dcp at 1he trailing edge (determined 

c 

graphically frcxn a faired plot of against cp) . The corre- 

sponding contoibution to the potential ;1b 

The velocity potential $ =^4 •’■‘5 o constitutes 

the net effect of putting the airfoil in the cascade (that is, the 
net interference effect) may now be determined by simple addition 
of the three ccmponents. Presumably, since tho calculations were 
made with Iho correct should bo the difference between 

and 

In the final step, is differentiated with respect to 

distance along the airfoil to get the corresponding Interference 
effect on the velocity v^ Tdilch should be the difference between 
Vj, and v^,. Convenient procedures for performing these calculations 
are discussed in the section entitled "Computational Methods." 
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ITERATION METHOD 


In the preceding sections the hasic ccaicepts and procedures 
of cascade interference calculations have "been outlined. In the 
present section, the application of such calculations in the pro- 
posed iteration method of solving cascade flow will he discussed. 

As first attempted, the method was essentially as follows: 

In the first step, is assumed to equal 0 ^ and a first 

approximation to is calciAlated on this hasis hy the methods 

just descrlhed. In the second step, is assumed equal to the 

sum of and this first approximation to and a second 

approximation to is computed. The succeeding steps follow 

the same pattern and are' continued until two successive 

& 

dlstrihutiona are essentially the same . The source-vortex method 
was used for the earlier approximations, hut the final approximation, 
when convergence is practically cojaplete, was made hy the contour- ' 
Integral method. This procedure, however, was found to converge 
relatively slowly in acme cases j and the general practicahility of 
the interference method depends on a sli^t modification of the 
source -vortex method. 


The modification depends upon "ttie ohservations that the 
contrihutlon of the eourcos and sinks to, changes hy relatively 

little from one approximation to the next and that the contrihution 
cxf the vortices to is nearly proportional to their total 

strength and relatively independent of their distrihution. Obviously, 
if it were exactly true that the contrihutlon cf the sources and 
sinks is constant and that the contrihution of the vortices is 
proportional to their total strength, only one interference calculation 
would he required and 15 ie solution could then he obtained ttoou^ a 
simple algebraic equation. Thus, let 

total circTilation on airfoil in. cascade 


8.C 


total circulation on isolated airfoil at same angle of attack - 
additional circulation 

constant contrihution of sources and sinks to T- 


Tg contrihution of vortices to Pg, when is assumed on all 

external airfoils 
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Then, "by the preceding as simp t ions. 


whence 


^T - J'f Tag 



T 


1 


( 1 ) 


Since the assimptions are not exactly true, the value of Prji 

so calculated is correapcoidlngly inexact; however, it is much 
closer to the true value than if 'it were taken simply as 
Pf + Pa Ta • Correspondingly, the potential 

S V 

+ 

^f + *ag + V* 


— ^a 
Vf V 


is much more accurate than the sum 


The second approximation is similarly adjusted. Thus, 
corresponding to. the 's&^-distrl'butlon just obtained, a new set of 

sources, sinks, and vortices are distributed along the mean line, 
and new values of P^^ and ^ay are calculated. Adjustment- 

follows, as before, frcm the equation 






where the subscripts 1 and 2 refer to the first and second 
approximations, respectiveiy , Solution for Pm gives 



^f + Ta 


■«v 



1 


( 2 ) 
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and, finally, tiie potential is given ty 






\ 


This simple modification of tlie procedure is so effective 
that in the cases tried, the first step gave solutions that would 
he satisfactory for many purposes and the procedure had practically 
converged at the second step. The additional complication of 
keeping the source-sink and the vortex effects separate so that 
and can he separately computed is relatively minor and 

amply repaid hy the rapidity of convergence. 

After the source -vortex method has essentially converged, a 
final approximation hy the contour -integral method is desirahle. 

In the cases computed, however, this final step was found to 
introduce only minor changes in the result. 


THE FLOW AT OTHER ANGLES OF ATTACK 


From a known velocity distribution at a given angle of attack, 
the angle of zero lift and the elope of the lift curve, together 
with the velocity distribution at any other angle of attack, may he 
obtained. For this purpose, the lattice is conveniently considered 
to he related conformally to an Isolated circle hy a periodic 
transformation, -vdiich mi^t he, say, of the type used in reference 6, 
8, or 9. The eaplicit form of the transfoimation, however, is not 
needed for the present purpose. 

The flow function in the circle (t ) plane that corresponds 
to the desired flow in the physical (z) plane is 


W = 


-—\e 
2n \ 


logg 


t + 
t - 


+ 


iXo 

e 


logg 



IT , 


( 3 ) 
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In the t -plane this flov may he interpreted, as that due to. the 
system of sources, sJnks, and vortices shojjn in figure 1. The 

■unit circle t = e^^ is a streaanline of 'the flov and "the circu- 
lation about any contour enclosing this circle hut not enclosing 

ihe points ^ = le^ is P (positive clockwise) . 


In the physical (z) plane, the complex velocities at 'the 
points z = <» and z are determined hy equation (3) and 

ihe transformation. Thus, 



and 


<iz 


=-v 


-iOQ 


- l^e . 
2d . 




-iOg 


where the angles and velocities aro defined in figure 1. The 
flow far from the la'ttico Is seen to he the sairg as -that of an 

-ICCq 

inf ini •be vortex row In ■the uniform flow ~^q® • should he 

noted (fig. l) that Xq = Oq + P, Xj^t=a|j^ + P, and Xg = ttg + p. 

In the following paragraphs ■ it will ho shown how to oh'fcain from ' 
the gl'ven spluticmi in cascade the parameter K and the stagnatiqaa 
‘points and 6.J. for the corresponding flow about- the. cirdo » 

These values fix the angle of zero lift and the slope of .the lif ■{> 
curve of the airfoil in cascade; together with the knoTO potential 
distribution 'bhey determine the conformal correspondence between 
the profile and the circle and, hence, -the velocity distribution 
at any angle of. attack . ; 

Since the airfoil contour (z -plane) is conformally related to 
the iinit circle (t-plane), it follows that at any given angle of 
at-tack' the change of .velocity potential from nose: to .tail 

stagnation point on bcth upper and lower sturf aces must he the same 
for -bhe circle and for the profile in cascade. These potential 
changes can readily ho obtained for the single solution on the 
lattice from the final ^^-distribution. The velocity potential. 

on the unit circle is. obtained from equation (3), Thus, 
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3^ 

2rta 


cos Xq l°8e 


cosh K + COB 0 


, cosli E - cos 0 


-V 


2 sin X, 


0 


tan 


'^ 8 In 0 
sinh K 


+ il- tan'^fe-0_ 

Yq&. tarJi K 


ih) 


and the change of potential frcm nose stagnation point 0^^ to tail 
stagnation point 0^ is 



1 

2jtcj 


<' cos Xq logg 

L- 


(cosh K - cos 0-t)(cosh K + cos 0 q) 
(cosh K + cos 0^) (cosh K - cos 0j^) 


+ 


2 sin Xq 


tan 


(sin 0^ - sin Q-^) sinh K 
sinh% + sin 0j^ sin 0^ _ 


+ 


VQd 


tan 


-1 


(tan 0y^ - tan 0 -j.) tonh K 
tanh^ + tan 0jj tan 0^ 


(5) 


This potential change may he obtained for either the upper or the 
lower surface. Two values are obtained depending on the choice of 
quadrant for the third term of equation ( 5 ). The condition of zero 
velocity at nose and tail stagnation points is 


sin 0 cos Xq - cos 0 tanh K sin Xq 


r 

27Qd 


sinh K = 0 


( 6 ) 


By Tise of the known values of P, and Xq, equations ( 5 ) 

and (6) can be solved simultaneously for 0j^, and E. 

Equation (6) can be considered as a quadratic in sin 0 and with 
an asstaned value of E determines corresponding values of 0j^ 
and 0^. Equation ( 5 ) then detemines By the proper choice 

of values of E, a curve of against E may be plotted such 

that at a point on this curve valiie of E at this 

■ point is the desired value ; the corresponding values of 0^^^ and 0^ 



ara tiwn given Tdj- eqmticm (6) . A convenient initial choice for E la the value that 
correBUonds to a lattice of eti'ai^t lines of the san© stagger and of about 10 percent or 
20 percent higher solidily. ligure 3 is of aid in this respect. The conqjuted values 
of K and 0^, together with equation (6), determine the angle of zero lift = 0) 

with respect to the airfoil chcrd, thtis, i 




-1 

tan 


tan 0-t Q 
tanh E ' ^ 


( 7 ) 


and the slope of the lift curve, hesed on mean velocity, is obtained hy differentiating 
equation (6) with respect to \qJ thus, 


dc^ !(. 

dn^ a Blnh E cosh E 


.( 8 ) 


A correspondence between points on the airfoil and points cn the unit circle may 
be obtained by comparing the values of ccaaputed by equation (it) with the values 

of 5>qi frcm the known potential distribution. The points (i,y) on the profile 

and 0 On the circle for which = 0^ are corresponding points. The velocity on the 

lattice profile for the Btream angle Xq< is ■ 


( -\ = 

ldt 

UaI 

dz- 


1 

vwvOT m m I JV I 


da 


idz 


co s cosh E (s in 0 - sin 0^) - sin Xq sinh E(cos 0 - cos 0^) 
- cosh%: - eos^0 ■ 


•vhere the teim in brackets, which represents the velocity on the, circle boundary, is 
obtained by differentiating eqliation (h) . It follows that the velocity corresponding 
to a new stream angle Xq* is 


H 


r 


HACA THf No. 1252 





coa cosh K( s3n 0 - Bin 0^) - sin alnli K(cos 0 - cob 0^) 
coa Xq cosh E(sin 0 ^ aln 0^) - Bin Xq Binh E(c6b 0 • cos 0^) 


( 9 ) 


The following relations, which ioacribe the flow far away frcm the lattice, are of interest. 
The Btream angles . X^ and. Xg at a = oo and a = -oo are 


Bln X^ 


+ JL 


\ 4*ow^ 

/V 

1 


>1 0 


COB X 


0 


aad 




X. = tan' 


cos Xr 


and the angle through which the fluid is turned by th'.- D.attlce is given by 


CD = tan 


-1 Tpd 


COB Xf 


- 

[S7^6j 


The rise in static pressirre across the lattice is 

An /Vi\2 /Vo'^2 




cos^Xq (sec^^ - aoc \) 


M 

\JI 
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EEMAEIKS ON CONTOUR MODIFICATIONS CORRESPONDING TO LOCAL 

PRESSURE CHANGES ■ 


In refersnce 14 , the modification of an airfoil contour to 
obtain, approximateli', desired email changes in the pressure 
distribution is discussed. The method, based on the formulas of 
reference 7, evaluates a sli^t modification of the confonaa^ 
transformation of the ci2rcle to the airfoil, such that the 
stretching factor at every point is changed in proportion to the 
desired relative change in local velocity. 

Although in reference l 4 the airfoil was assumed to lie in 
a straight uniform field, the treatment ls_equally applicable 'vdien 
the airfoil is in a curved or distorted flow field. Accordingly, 
the procedxire should be applicable to airfoils in cascade, 
provided the same modification of the orternal airfoils leaves 
the disturbance flow field essentially unaffected. This condition 
may not always be satisfied; however*, in such cases the method 
could possibly be improved by a procedure analogous to that 
described in the section of the present paper entitled "Iteration 
Method." ■ 


COMPUTATIONAL METHODS 


The basic theory has been presented. In the following sectioiB 
some oi’ -the methods used for performing the actual computations will 
be discussed. 

Sele ction of points for evaluation of distttrbance flow . - The 
determination of the oompenBatlns flow by the method of reference 13 
requires that the disturbance flow be evaluated at points "uhat, by 
the conformal transformation, correspond to points equally spaced 
about the circle. These points, which are located by reference to 
the confoi’mal transformation, are preferably chosen so that one is 
at the trailing edge. Experience has shown that, for the 
prelimiriBry approximations, 12 points at 30° intervals yield 
acceptable results. In the final step by tho contour -integral method 
the use of 24 points is preferable in order to Improve accuracy, 
especially near the nose. An acceptable compromise is to evaluate 
and directly for only tiie additional points that are near 

the leading edge and to pick off the values at the other additional 
points from a faired curve . 
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Inasmuch as Tralues for the 12 -point and 24 -point methode are 
not included in reference 13 the following table is presented: 


k 

c 



n = 6 

n = 12 

1 

0.62201 

0.63298 

3 

.16667 

.20118 

5 

.04466 

.10860 

T 


. .06394 

9 

— 

.03452 

11 


.01097 

. 


Evaluation of and Integration of eqmtion (36) of 

reference T along tiie circle boundary yields the values of the 
potential at points on the airfoil ^ follows: 

^f ^ ^0 
— == 2ae 


J 


<p Bln(ci + g) - cos(a + cp) 


( 10 ) 


angle of attack 
angle of attack for zero lift 

radius of the circle to which the airfoil transforms 

angular position along the’ circle’, as detenained by the 
tr ensf ormati on 

If the transformation has- beoh performed as recommended in 
reference T, the ccsnstant (a) will be slightly less than one' fourth 
the chord. Although the ' potential discontinuity (coi’respond.ing to 
the circulation) may, without loss of generality, be placed at any 
point on the contour, the trailing edge will generally be found to 
be the most convenient location. ■ ■■. 


where 

a 

3 


ao 

cp 


^^0 
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The additional Telocity t„ 16 given "by the deri^fative along 
the 'sulrface ■ it may he determined hy graphically differen- 
tiating with respect to the circle angle cp and multiplying 

this Slone hy Thus. 

ds . 

Tq * ,Vo<ia ’ V^dq) ds ' (ll) 


The value of 
reference 7* 


^ may he obtained from equations (37) aJid (38) of 
Thus^ ■ ■ 



or 



( 12 ) 


■where the symhols e, 6, and, are defined in reference 7- 

The cascade solidity need he taken in-to account only when the 
airfoil sketch to he used with figure 2 .is_ /constructed; For the 
Buhse'quent c.eJLc^1lations, any convenient airfoil, chord may he used, 
provided only. -fchat the- same chord is used for the ex-temal airfoils 
and for ■fche cen'tral airfoil.- The reason is as follows; The 
strengths of the aingulari ties used to represent the external ■ 
airi’olls are proportionai. to the assumed airfoil chord; hence the ' 
additional potentials induced on the central airfoil will be 
proportional to the assumed, chord. Since both the additional 
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potential and the distance s ^eiosig, the contour s^e pro- 

portional to the chord, the additional velocfrby v„ = -3#- will 

Gr 


ds 


he independent of the chord. 


The chord may then conveniently he chosen as that corresponding 
to a veO-iie a = 1 since a wovild then not appear in equations (lO) 
and (12) 


The net velocity at a. point on the airfoil surface is the 
algebraic sum of the velocity on the isolated airfoil and the 
induced velocity Vg^ at that point. 


Selection of voi-ticea for source-vortex method. - For cascades 
of about unit solidity, the vortex distribution for an airfoil of 
conventional design may be represented by five vortices spaced on 
the mean line at O.l, 0 . 3 ^ 0 -/, 0 . 7 , and 0.9 of the chord. The 
sti^ngths of the vortices are detemlnod by liie known chordwise 
distribution of potential on the upper and lower surfaces for 


the given approi.Lmation . The difference in potential between the 
upper and lower surfaces at 0.2 chord is thus approximately the 
total vorticity between the le-ad?.ng edge and 0.2 chord and Is 
considered to be concentrated in the vortex at 0.1 chord; similarly, 
the increase in this potential, difference between 0.2 chord and 
0.i^ chord yields the strength of the vortex at 0.3 chord, and so on. 


The total vortex strength miist satisfy the equation 




S election of sources and sinks for source-vortex method . - The 
selection of sources and sinks to represent the thickness distribu- 
tion of aiifoils is less readily systematized than is the selection 
of vortices to represent the lift dietributicfn. For conventional 
airfoils, a reasonably satisfactory representation is geherally 
attainable with a source at about 0.025 chord, a sec end source 
midway between the nose and the position of max.imum thickness, and 
sinks at 0.5, 0 . 7 , and 0.9 of tlie chord. The strength of each 
source or sink Is taken as the difference between the "inteirnal flow" 
at a station midway between it and the preceding source, and the 
internal flow at a station midway betireen it and the following 
soTirce . This internal flow at a given station is estimated to be 
the product of the thickness and the average of the upper and lower 
surface velocities at that station. 



20 


NACA TN No. 1252 


OlDTiotxaly, not all airfoil shapes vill he best treated 
according to the pattern just described; however, little ingenuity 
is required to adjust the treataont to a parti cular shape . In any 
case, the total source strength must equal the total sink strength. 


PROCEDURE 


A suggested step-by-step procedure is as follows: 

(1) Obtain the velocities on the airfoil at ihe given angle 

of attack in a uniform stream by the method of reforence.T* This 
step also determines a conformal correspondence between points 
(r, y) on the airfoil and angles cp on a circle, and hence the 
potential distribution by equation (lO) . 

( 2 ) Using the procedure described in the section entitled 
"Ccmiputational Methods" choose sources, sinks, and vortices to 
represent the airfoil. 

( 3 ) Choose points around the airfoil at which the distur- 
bance function is to be found; these points aro conveniently 

chosen, by reference to the conformal transformation, to correspond 
to 12 equal intervals about the circle. By use of figure 2, 
determine at these points the contributions, to and of 

each source and vortex row. Sum separately the values due to 
sources and vortices at each point. 


( 4 ) Form the ccanpensating functions both for 

vortices and sources and determine the conjugate functions by the 


method of reference 13 . Plot 


slope at the tarailing-edge point, 
determines the circulation changes 


against , p and measure 

The -relation P = -2« / 

\ / T .E . 


and 




due to the 


source and vortex rows. Obtain P^ by means of equation (l) . 


T 

( 5 ) At each point 

(a) Sum the values of Oa. 


and Or 


duo to the 
Pm 

vortex rows and multiply by the ratio 

Pf 

(b) Sum the values of 0^^ and due to the rows 

cf sources and sinks. 

(c) Find ^ - 
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( 6 ) Sum the terms (a), (h), end (c) of step ( 5 ) to get $ 5 
plot against -ttie circle angle cp, and msasui’e the elopes at 

the points used In the original conformal transformation (step (l)) 

at which noints the stretching factor will he known. The 

additional velocity is given by equation (ll); the not velocity on 
the airfoil surface is the sum. of the additional velocity and the 
velocity on the isolated airfoil. The corresponding total 
potential is 5,^, = > where is known from 

step (? ) . 


Using this new potential, and velocity distribution, repeat 
the procedure, starting from step (2) . The only modification is 
that r,p (step (4)) is now obtained from equation ( 2 ), anrl in. 

step ( 5 a) the correction factor is J^t^/^Ti* Th® process is 

continued imtil the changes in lift and velocity distribution 
become small. For practical purposes, the results obtained in this 
manner may be entirely satisfactory. More acdurate results may be 
obtained, however, by application of the contour-integral method 
as described in the following three steps , 


(7) Locate the points on the airfoil that correspond, by the 
conformal transfonna.tlon, to points midway between those already 
located in step (3) • Place the airfoil drawing on figure 2 with 
the origin, in turn, at each of the 12 points at which values, are 
known from step ( 6 ) (considered as z-points)^ and read the chart 
at each of the 2k points (considered as z’ -points). As previoixsly 
noted, seme of these points may be neglected. For each of the 24 
(or fewer) points plot the 12 values of ^ read at that point 
against the .12 corresponding values of By planimetry find 

the area between the faired curve and the , ^^-ajcis to determine 5 ^^. 

The value of is determined eimi.larly frem a plot of the 12 

values of against the coiresponding values of 


( 8 ) Form the function determine its ccaijugate 


the circulation change is 


c. 


^■n - r_ 


cp 

7 ^' 


-2rt[ 


and the potential 


(9) Sum the terms and to get plot against 

the circle angle cp, and measure the slopes. The velocities on 
the airfoil surface in cascade are obtained as described in step ( 6 ) 
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Unless this volocity dlstrlMition differs widely from that obtained 
in the preceding anproxiitiation. It should not be necessary to 
repeat the procedure 

The velocity distribution at another angle of attack may be 
obtained as follows: 

(a) Solve equations (5) and (6) for 9^^, 6^, and K. A 
method of solution is Indicated in the discussion following 
equation (6). The angle of zero lift and slope of the lift, curve 
may then be obtained from equaticaiB (7) and (8). 

(b) Obtain the potential distribution ^ as a function 

of G (equation compare with the known $qi to get a 

correspondence between 0 and position on the airfoil. Equation (9) 
then yields the velocity distribution at stream angle . 


niUBTRATIVE EXAMPLES 


Example 1 . - The velocity distribution was obtained on the 
MCA i(-4l2 airzfoil in the conf igua’ation shown in figure k, where 
p = 0°, a = 1.032, and Xq = 9.7°- This example has been 

treated in reference 8. In accordance with the foregoing procedure, 
results as follows were obtained: 


(1) In figure 5 is shown the chordwise velocity distributions 
of the isolated, airfoil at the angle of attack of 9-7°^ as obtained 
in a second approximation by the method of reference 7- The lift 


coefficient at this angl.e of attack is I.67 (that is, 




0.837), 


the angle of zero lift of the airfoil is -4.24°, eind the slope of 
the lift curve is 6.95 radian. 


(2) By use of the procedure suggested in_ the section entitled 
"Ccaaputatlonal Methods, five vortices, two sources, and three sinks 
were chosen to represent the airfoil Initially (fig. 6 and table I). 


(3) With the first location at the trailing edge, 12 locations 
on the airfoil were found corresponding to 30° intervals of the 
ci3:cle angle cp. These locations are shown in figure 6. (The 
primed points correspond to 15° intervals!) Eeadings taken at 
these points from figure 2 are given in. table II. TliQse readings, 
multiplied by the appropriate source a^ vortex str^gths, yielded 
the values of and due to sources end vortices given in 

table III, 
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(i}-) The conjtigate fimctlons 0^, were determined hy the 12- 
I>oint method, and are given in table 17. The slopes of these 
functions at the trailing-edgo point yielded circulation changes 

= 0.006 and ^ = -O.53S, from which (eauation (l)) 

C /r 


cTo 


0 


■=^= 0.513. This value corresponds to a first approximate lift 
°'^0 


coefficient in cascade (c-^ = I.03) . 

( 5 ) In table 17 are given the values of 


'^d-tr 


vortex rows multiplied by the ratio — - (equation (l)), the values 

' ^ f 

of and §4 due to source rows, and the function 


$p = (Tt ■ Pf ) 2rt 


(6) The additional potential 0a=’0j_ + ^o+0n is plotted 

in figure 7. Slopes of this function were measured at points at 
which the stretching factor is known from step (l) . The additional 
velocity Vg^ was -then computed bj eqmtion (l0)j the algebraic 

sum of Vg and the velocity in isolated flow yielded the cascade 

velocity (fig. 5) . This velocity distribution, together with the 
total potential formed the basis for a second anproxlmation 

p 

(figs. 5 and 7). Results of this approximation are = O.OO6, 

cVo 


fay 

- - 0 . 365 , and c^ = 0.99. CcaiparlBon of the velocity 

distribution with that of the first approximation shows that the 
process has satisfactorily converged. 


(7) The same 12 points around the airfoil were chosen as 
z-points; these, togot.her with four others at 15° intervals around 
the nose (primed points in fig. 6) were used as z’ -points. 

Readings from the chart (fig. 2) are given in table 7. These values 
were plotted against total potential j (arbitrarily fixed at 0 
on the lower surface- at the trailing edge). (A sample curve is shown 
in fig. 8.) These curves were integrated by plemlmetry. The 
resTolts - the disturbance potentied-s and stream functions 0 ^ 

and - are given in table 71, 
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(8) Tfae fimction (table Yl) was obtained by 24 -point 

0 

harmonic analysis and synthesis, \7ith the use of interpolated 
values of for the points at which it was not found explicitly. 

The slope of the curve at the trailing-edge point yielded 

= :-0.34^+, from ■vdiich a lift coefficient c, = 0.99 was 
cVo ^ ^ ^ 

obtained. 

(9) The additional potential + §p is plotted 

in figure 7- The velocity distribution was obtained as before and 
is plotted in figure 5* The ‘process appeal's to have essentially 
converged. . _ ... _ 


Simul'baneous solution of eqjip.tions (5) and (6) (table VII) -fco . 


find the value of K at which 


A0„ 


cVq‘ ".cV, 


0 


gave K = 0 . 3033 , 


«n = 


*7.?7°, and 0^ = 181 . 72 ^. Sguatians ( 7 ) and (8) -then 


yielded the angle of zero lift = -5-75° 'tlia slope of the 

dc. 


lift curve = 3 . 71 . 

aaQ 


These values may be ccmpared with 


q = - 0 . 94 ° and 


dc 


= 3*71 from reference 8. 


In figure 9 is shown a plot of -the potential against 9 

computed by equation (4) . A constant has been added to make "the 
potential eq’jial -bo zero on the lower siurfece at the -trailing edge. 
The Icnown total po-fcential in cascade and -the corresponding 

values of x/c are given in -table VTII. Values of 0, picked off 
the plot at points -ifhere is equal to the given -values of. 

are shown in the adjacent column. The correspondence between 
airfoil posltl.on and the angle 9 is thus determined. For the 
flow angles 1 q* = 1.8l° and Xq’ = -3-9k^, the velocity 

distributions were caupu-tsd by equation (9). In figure 10 these 
results are compared wi-th -the distributions given in reference 8. 
The main results of the calculations are summarized in table IX. 


E xample I I ■ - In an effort to ob -tain in the slmplost possible 
manner.a reference solution at lai*^ bla^ angle, conoeming the 
accui’acy of which there could be little doubt, a lattice was 
derived by a modified Joukowski transformation, Thi.s trans- 
formation Is discussed in' detail in the appendix. The cascade 
configuration is sho-wn in figure 11 where p = 45°, o I. 006 , 
and Xo = 49°. This lattice will be referred -bo as the "derived 

M 

airfoil lattice . 
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The procedurs followed for -She eom’ce -vortex method was slmller 
to that of the first example: the caJ culat5.one are outlin<id in 
figures 32 to li^. Becanss of 'I'he unusual shapo of this profile, 
only one source was used and an additional sink was inserted at O .3 
chord (fig. 12). Frau a lift coefficient = 0.8^4- in isolated 

flow, a single approximation yielded a lift coefficient c^^ = 0.5^ 

in cascade, which was the same as that derived from the solution hy 
conformal ti'ansfoimation. Since the computed changes in vortex 
distribution were sma3.1, no further approximations were made by 
this method. By reference to the velocity distribution of this 
approximation (fig. 13 ) ^ the process may be seen to have ess<^ntleUy 
converged to the correct solution. 

The final contoiir integration resulted in a lift coefficient 
c^ = 0 . 5 ^ and the velocity distribution shown in figure 13- The 
main restuLte of the calculations are sirnmarized in table X. 


Lang3.ey Ifemorlal Aeronautical Laboratory 

National Advisoiy Ccmmittee for Aeronautics 
Langley Field, Va., January 1C, 19^7 
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APPENDIX 

DERITiilD AIRFOIL LATTICE 


The syrnhols used in the appendix are defined in figure I 5 
and should not he confused with similar symbols used in the main 
text of the paper. 


Consider the transf o 2 miation (reference 10), 


z = 


-13 


2it\f - eK. 


t + eE IB . t + e-E 

+ © ^ logg 


t 


- e-E, 


(Al) 


The unit circle (^-plane) heccaneB a lattice of horizontal 
straight lines in the z -plane, spaced at unit Intervals along the 

stagger line , making an angle ^ - 3 wJ th the axis .of reals . The 

solidity of this lattice is 


a » f cos 3 log. 


'/ ainh^E + cos^P + cos 3 
slnh E 


+ sin p tan 


-1 


-g-ln .3 


I'^inh^ + coB^ 


This relation is plotted in flgiire 3* 

A closed curve enclosing the points ^ = te"^ hut not 
enclosing the points ^ will transform hy equation (Al) into 

an infinite lattice of closed shapes in the z-plane, spaced in the same 
manner as the strai^t-line lattice. Such a curve is the circle 

iS 

= e + re 


= l.OTe 


id 

^ + 0 . 


09e 


in 

3.75 


t 
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This circle, where (3 = ^5° and K = 0.331, ’becomes -the lattice of 
profiles that has "been referred to as the derived airfoil lattice. 

A flow for which this circle is a strecmline and which, in the 
z-plane, has no singularities outside the profiles, is that due to' 
the system of sources, sinks, and vortices shown in figure I 5 . The 
velocity on the circle "boimdary due to this system is 



= A cos Xq + B sin Xq 



where 


A s= e 




sin(^ - 61) Bin(i^ 

- cos(^ - Hg - cos(5^ - 


B = e 




.£ 1 . 


H. 


- cqb((^ - 6^) Hg - cosC^if - 6g) 


C = e 


-Vo 




- cosC^f - O^.) ^2 - cos(^ - Sg) 


and 


5^ = tan 


-1 r sin 5 


e^ - r COB 6 


5 = tan 

2 


-1 -r sin 6 


+ r COB 6 


H = Ifm + i 


J = - m'' 


m^ = e - 2re^ cos 6 mg = e \/r^ + + 2re^ cos 5 


5he constant P which is tho circulation about each profile 
(positive clockwise), is determined by the trailing-edge condition as 


r 

T 




cos Xq -f- — sin X, 


1 


(A2) 
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'where A, B, and C are evaluated at the angle ^ -which 
correspcmds to -fche trailing edge of the -profile . The angle of 
zero lift t) -with respect to -fche airfoil chord, is oh-fcalned 
frcm equation (A2) by setting P = 0; -thus, 

- 1 . 

11 = -tan ^ - 3 

The stretching factor from the circle to the lattice is 



where 


D =1 [^cosh 2K - cosh 2"4/ cos 29)^ + (sinh 2^ sin 20)^ | 
^ cos^p cosh^K(cosh^ - cos^0) 


1 


+ 1+ sin^p slnh^Ccosh^-^ - sin^9) - sin 2p sin 20 sinh aCj 

and •({/ and 0 are obtained from ij/Q, r, and 8 as 

^0 

Q _ tan"^ e sin 0 + r sin 5 
e ^ cos 9 + r cos 8 


e^ = r cos(0 - 5) 


. v4^o . 


sin ^ (0 - 5) 


The velocity at any point on the surface of a profile is 
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STREHGTHS OF SODRCffiS /HD ■VORTICES CHDSHN TO BEPKESBRT 
THE BfACA 4432 A3HFOIL LATTICE 



"^0 

Source 
location 
(fie* 6) 

m 

Vortex location 
(fig. 6) 

First approii- 
matloQ 

Second approxL- 
oatlcnx 

First approxi- 
mation 

Second approxi- 
mation 

P 

0.379 

0.290 

a 

0.097 

0.101 

6 

.184 

.098 

r 

.044 

.044 


lOP 

r\£/^ 

£ 

aL a 
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-.045 

-.047 

1 
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"•053 
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TABLE n 


CHART EBADIHGS FOR FACA^OS AIRFOIL LATTICS, SCUmCB- (S3 

YOHTKI MSmOD 



9 for vort«i rov of unit 
strength 


f 

for Tortei row of mj 
BtMngth 

t 


^\tolgjn at 
\(flg- 6)- 

Headlng 
(fig. 6)“ 

^7 

a 

r 

P 

Aff 

7 

6 

t 

iO 

- Cl - ii ' 

t 

4 

ct 

P 

7 


c 

i 

n 

a 

O.OCE 

0.006 

0.008 

ou}n 

0.010 

0.004 

0.001 

-0.2U 

-0.184 

-0.176 

-OJ20 

-OMTO 

-0.025 

-0.003 

b 

.002 

.008 

.008 

.011 

.030 

.003 

0 

-.182 

-.158 

-.150 

-.096 

-.04? 

-.013 

0 

c 

.OOli 

.008 

.008 

.009 

•005 

-.001 

-.002 

-.115 

-.092 

-.083 

-.044 

-.011 

0 

-.010 

d 

.005 

.006 

.006 

.004 

-.002 

-.007 

-.007 

-.045 

-.030 

-.026 

-.004 

0 

-.019 

-.056 

e 

.002 

.002 

.002 

-.004 

-.012 

-.013 

-.013 

-.OCrf 

-.001 

0. 

-.002 

-.025 

-.040 

-.125 

f 

0 

-.001 

-.002 

-.006 

-.015 

-.033 

-.011 

0 

-.001 

-.001 

-.019 

-.057 

-.U5 

-.174 

8 

0 

0 

0 

-.003 

-.005 

-.002 

0 

0 

-.003 

-.002 

-.024 

-.063 

-.120 

-.183 

h 

-.002 

0 

.001 

.002 

.005 

.010 

.016 

0 

.001 

0 

-.012 

-.044 

-.090 

-.156 

1 

-.012 

-.007 

-.006 

.001 

.008 

.015 

.023 

-.014 

-.006 

-.004 

.001 

-.013 

-.048 

-.096 

J 

-.oeo 

-.015 

-.ow 

-.006 

0 

.010 

.016 

-.O5B 

-.{«2 

-.037 

-•010 

•001 

“ *010 
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fc 

-.016 

-.011 

-.010 


.002 

0 

.004 

-.125 

-.102 

-.094 

-.050 

-.018 

0 

-.006 

1 

-.005 

0 

0 

.005 

.004 


0 

-.105 

-.161 

'.153 

-.097 

-.048 

-.015 
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TABLE m 


CCHOBIBTJTICRS OF IHDIVIIJUAL SODECE AHD VDHITEI BOWS TO THE 
DISTDBBAHCE FLOW FOTCTiaf OS TBS HACA AIEFOIL 3 H 
CASCALEj FIEST APEROrEMATICir, S0DECS-V(SiTK3C MBiaOD 



P 

8 

e 

£ 

n 

a 

7 

e 

— 

£ 

4 











Beetdlng 

at- 


due to Tortex rows 

dais 

to source rovs 


( fig . 6 )- 











a 

0.0030 

0.0020 

0.0013 

0.0004 

0 

0.0205 

0.0077 

-0.0030 

-0.0011 

- 0.0001 

t 

.0030 

.0020 

.0013 

.0003 

0 

.0176 

.0066 

-.0020 

-.0006 

0 

c 

.0030 

.0016 

.0006 

-.0001 

-.0001 

. 0 U 2 

.0036 

-.0005 

0 

-.0005 

d 

.0023 

.0007 

-.0002 

-.0007 

-.0004 

.0044 

.0011 

0 

-.0008 

-.0030 

e 

.0008 

-.6007 

-.0015 

-.0013 

-.0007 

.0007 

0 

-.0011 

-.0018 

-. 00 ^ 

f 

-.OOCk 

-.0015 

-.0019 

-.0013 

-.0006 

0 

.0001 

-.0024 

-.0052 

-,oaS^ 

g 

0 

-.0006 

-.0006J 

-.0002 

0 

0 

.0001 

-.0027 

-.00^ 

-.0097 

h 

0 

.0005 

.0006 

.0010 

.0008 

0 

0 

-.0019 

-.0040 

-.0083 

1 

/-.0026 

.0002 

.0010 

.0014 

.0012 

.0014 

.0002 

-.0006 

-.0022 

-.0052 

J 

-.0057 

-.0011 

0 

.0010 

.0008 

.0056 

.0016 

0 

-.0004 

-.0022 

k 

-.0042 

-.0009 

.0003 

0 

.0002 

.0121 

.0041 

-.0008 

O' 

-.0003 

1 

0 

.0009 

.0005 

.0001 

0 



.0180 

.0067 

-.0021 

-.0007 

0 


due to TortOi rows 

due 

to source rows 


a 

-0.0698 

- 0.0221 

- 0.0090 

- 0.0024 

- 0.0001 

0.0002 

0.0004 

- 0.0004 

- 6.0002 

0 

1) 

-.0599 

- .0176 

-.0060 

-.0013 

0 

.0002 

.0004 

-.0004 

-.0001 

9 

c 

-.0348 

-.0081 

-.0014 

0 

-.0005 

.0004 

.0004 

-.0002 

6 

.0001 

d 

-. 0 U 4 

-.0007 

0 

-.0018 

-.0029 

.0005 

.0003 

.0001 

.0003 

.0004 

e 

-.0005 

-.0005 

-.0032 

-.0039 

-. 00 ^ 

.0002 

.0001 

.0005 

.0006 

.0007 

f 

-.0003 

-.0035 

-.0073 

-.0111 

.-.0091 

0 

.0001 

.0006 

.0006 

.0006 

g 

-.0010 

-.0044 

-.0081 

-.0116 

-. 00 ^ 

0 

0 

.0002 

.0001 

0 

h 

.0002 

-.0021 

-.0056 

-.0087 

-.0081 

-.0002 

0 

-.0002 

-.0004 

-.0008 

1 

-.0023 

-.0002 

-.0017 

-.0046 

-.0051 

-.0012 

-.0003 

-.0003 

-.0007 

-.0012 

i 

-.1057 

-.0018 

.0001 

-.0009 

-.0022 

-.0019 

-.0006 

0 

-.0004 

-.0008 

k 

-.0386 

-.0092 

-.0022 

0 

-.0003 

-.0016 

-.0004 

-.0001 

6 

-.0002 

1 

-.0610 

-.0178 

-.0061 

-.0014 

0 

-.0005 

0 

-.0002 

6 

e 
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(i ... • ^ 

is ^ C 


Point 

"^0 


SonrcsB 

a o-dooo 0.001t 

ij jxxji .ooU 

c .0007 J )015 

a .0016 =0009 

0 .ooei - .0003 

f .0019 - .ooeo 

e .0003 - .0030 

ii -.0016 -.0032 

1 -.0037 - . 0 Q 3 P 

J -.0037 .0007 

-.0023 

-.0007 
.ooool .00 


TOiH. HnracT OF sooBCE am? vcEm M«B, asii coERESpcaramG 
AID TELOCITIBS, OK SAGA ATHTOIL Bt CASCAlBi FIHST 
SOCSCK^VCItTElX MEIISOD 

7olue& appoints aa^flgare 6 
^ ^ ^ * 7 . 

% *0 % (^av*' ^ 5a_ 


DERlViUJ POnBirnALS 
APERCraiATICO. 


. Vortices 

. ■ I , — 

0.0240 -0.1034 0.0013 D.ooer 
.0216 -.0848 .0414 .0066 
J3130 -.0448 .0513 .0050 
.0017 - .0168 .0314 ,0017 
-.0087 -.0146 .0040 -.0034 
-.0167 -.0313 -.0050 -.0057 

-.0177 -.0347 .0037 -.0014 
-.0142 -.CG43 .0044 .0029 

-.p64 -.0139 -UX)96 -00p2 
.0046 -.0205 - .0356 -ux>30 
.0151 -.0503 -.0505 
,(fil9 -.0663 -.0372 
.^40 -.1034 .0013 


- ait 

0JJ049 -0.0270 0.0033 


-.0065 

.0014 

.0045 

-.0053 

-.0236 

-.0338 

-.0218 


7 


CTq 

CTq 

'it '■ 

0.0033 

0,6056 D 

-.0015 

.5355 

-.0310 

.3592 

-.0646 

.1429 

-.1436 

-41365 - 

-.1872 

-.1I26 - 

-.2084 

-.0459 - 

-.2290 

.1642 - 

-41566 

Jt802 

-.2884 

.8365 

-.3136 

1.1558 

-.3218 

1.3716 1- 

-.3206 

14^51 L 


0.6089 

.5340 

•3282 

.0581 

-.1801 

-.2996 

-.2543 

-.0648 

.2236 

*5481: 


Talosa at points iihsre ^ le Imown 

u 5 

I ZS: h. 

c oVq dip 4 s 7g Tq 

I (a) [ (a) 

tfpper BiBcfaoe 


04)125 

-0.0356 7.153 

-0.254 

2487 

■ 2.033 

^fWiTtn 

- t ^kn 

.. lOh 

2 ., 0 ^ 

1 .000 

.1000 

-.0456 3.500 

-.160 

1.853 

1.693 

.2000 

-.0542 2.733 

-.148 

1.719 

1.571 

.4000 

-.0606 2418 

• 

0 

1.523 

1.383 

.6000 

-.0506 2.133 

-.108 

1.345 

1437 

.8000 

-.0226 2.496 

-.057 

1.178 

1.121 

.9000 

-.0078 3426 

-.025 

14)78 

1.053 


Lower elirfatfe 

0.0125 - 0.0451 8 j 263 - 'O .372 0 A 58 

.0500 -. 06514.603 -.300 -.314 

.1000 -.0847 3.304 -.200 -.541 

.2000 -.1027 S. 413 -.248 -.676 

■ -.221 



TwXocltiss a X m.ig the muXace ox'o 

Itao \ippar snrfBce at tha laTill-ln{; ed^. 


W MIT^ IT ATnrTQ/'tDV 
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TABLE V 

CHART HEADINGS FOR IHTSGBATICH WITH RESPECT TO 
HACA AIRFOIL IN CASCADE; 
CCBTODS-INTECRAL HEISOD 



a 

b 

c 

d. 

a 

f 

8 

h 

1 

j 

k 

1 


*1.08i8 


0.3106 

0.0401 

-0.2022 

-0.3280 

-0.2911 

-0.1064 

0.1783 

0.4986 

0.7954 

1.0060 

^0.58T5 



i . g ■ 1 

a 

0 

0 

•0 .001 

-0.0040 

-0.0095 

-0.008 

0.004 

0.020 

0.029 

0.0215 

0.006 

0.001 

1) 

0 

0 

0 

-.003 

-.0085 

-.006 

.005 

.021> 

.0275 

.020 

.0065 

0 

c 

-.0005 

0 

0 

-.0015 

-.006 

-.005 

.006 

.0195 

.0215 

.0105 

-.0005 

-.003 

d 

-.00^5 

-.003 

-.0015 

0 

-.0015 

-.001 

.0065 

.0145 

.010 

-.0035 

-.0125 

-.009 

6 

-.010 

-.008 

-.006 

-.0015 

0 

.0005 

.004 

.0055 

*.006 

-.0205 

-.023 

-.016 

e' 

-.011 

-.010 

-.007- 

-.002 

0 

0 

.0025 

0 

-.012 

-.025 

-.026 

-.018 

f 

-.009 

-.008 

-.005 

-.001 

.0005 

0 

.0005 

-.0025 

-.014 

-.026 

-.025 

-.014 

f' 

-.002 

-.001 

0 

.002 

.002 

0 

0 

-.0035 

-.014 

-.024 

-.0215 

-.010 

e 

.007 

.007 

.0065 

.0065 

.0045 

0 

0 

-.002 

-.010 

-.018 

-.015 

-.002 

b' 

.015 

.015 

.012 

.011 

.0055 

0 

0 

-.0005 

-.0065 

-.0125 

-.006 

.006 

h 

.OZk 

.022 

.019 

.0145 

.0050 

-.0035 

-.003 

0 

-.0025 

-.006 

.0005 

.015 

h' 

.028 


.022 

.014 

.001 

-.008 

-.0065 

.001 

-.0005 

-.0015 

.0055 

.0195 

1 

.0295 

.0285 

.022 

.0095 

-.006 

-.0155 

-.0115 

.003 

0 

.0005 

.006 

.0215 

J . 

.022 

.020 

.0105 

.0050 

-.0205 

-.027 

-.OIS5 

.006 

0 

0 

.0035 

.016 

k 

.008 

.0065 

-.0015 

-.012 

-.023 

-.025 

-.014 

-.0005 

.006 

.0035 

0 

.004 

1 

.001 

0 

-.003 

-.0085 

-.015 

-.017 

-.001 

-.015 

.021 

.0155 

.004 

0 


$ 

a 

0 

-0.002 

-0.024 

-0.083 

-0.154 

-0.205 

-0.217 

-0.186 

-0.130 

r0.066 

.0.018 

-0.002 

b 

-.002 

0 

-.012 

-.063 

-.132 

-.178 

-.187 

-.159 

-.102 

-.045 

-.0075 

0 

c 

-.023 

-.012 

0 

-.Cfi05 

-073 

-.111 

-.120 

-.095 

-.048 

-.0085 

.001 

-.0125 

d. 

-.061 

-.060 

-.0205 

0 

-.018 

-.042 

-.048 

-.030 

-.004 

.0025 

-.023 

-.062 

a 

-.153 

-.128 

-.070 

-.0175 

0 

-.006 

-.0075 

0 

.002 

-.023 

-.077 

-.131 

e' 

-.iSlt 

-157 

-.091 

-.031 

-.003 

-.0005 

-.0015 

.002 

-.0035 

-.040 

-.100 

-.159 

f 

-.204 

-.iSo 

-.109 

-.042 

-.007 

0 

0 

.0015 

-.ou 

-.053 

-.119 

-.181 

f' 

-.214 

-.189 

-.119 

-.048 

-.0065 

0 

0 

-.001 

-.017 

-.064 

-.129 

-.191 

8 

-.212 

-.186 

-.116 

-.047 

-.008 

0 

0 

-.0015 

-.0175 

-.063 

-.128 

-.189 

s' 

-.205 

-.176 

-.110 

-.0415 

-.005 

.0015 

0 

-.001 

-.0145 

-.056 

-.120 

-.181 

n 

-.184 

-.158 

-.093 

-.030 

0' 

.0015 

-.0015 

0 


-.044 

-.103 

-.160 

■9 

-.158 

-•133 

-.073 

-.017 

.0035 

-.002 

-.006 

-.002 


-.0295 

-.ote 

-.136 


-.127 

-.101 

-.048 

-.0035 

.0025 

-.011 

-.017 

-.0085 

0 

-.015 

-.058 

-.105 


-.062 

-.044 

-.008 

.003 

-.0225 

-.054 

-.063 

-.044 

-.0155 

0 . 

-.016 

-.047 


1 

b 


.001 

-.012 

-.081 

-.123 

-.130 

-.103 

-.059 

-.016 

0 

-.0085 

■ 

-.002 

0 

-.013 

-.009 

1 

H 

-.079 

-.188 

-.16^ 

-.108 

-.0475 

-.009 

0 


"tJpT)er eurfaoo at trallliig edge. 
- ^^Lover surface at trailing edge . 


H ATICB AL ADTISOEr 
CCMMITEKE FCE AIECHAOTICS 





36 


NACA TN No. 1252 


TABIX VI 

SSBXTED FOTEREEAIS AH) TBLOOZmS OH HACA 1(412 AIKTOIL II OASCADE; 
COITCOB^ilXGStAL MBIBOD 


Values at points on figure 6 

Values vbere 

g Is known 


B 

m 


m 

*P 

«a 

*f 

m 

X 

BE 

odq> 

a 








°^0 

“^0 


0 

oVq 

d 8 

wm 

a 


■ 











lol 



B 

0,0281 

-0.0687 

0.0021 

- 0.0144 

0.0158 

0.5717 

0.5875 

T^?per surface 

fl 

.0259 

-.0571 

.0271 

-.0432 

.0098 

.5302 

.5400 

0.0125 

- 0.0403 

7.153 

-0.288 

2.287 

1.999 

1 

.0164 

-.0314 

.0382. 

-.0720 

-.0174 

.3280 

.3106 

.0500 

.1000 

-.0496 

-.0525 

4.541 

3.500 

-.225 

-.184 

2.002 

1.853 

1.777 

1.669 

II 

.0033 

-.0098 

.0257 

-.1008 

-.0718 

.1119 

.0401 

.2000 

-.0667 

2.733 

-.182 

1.719 

1.537 

6 

-.0129 

-.0049 

.0090 

-.1296 

-.1335 

-.0687 

-.2022 

.4000 

.6000 

-.0650 

-.0452 

2.218 

2.133 

-.144 

-.096 

1.523 

1.345 

1.379 

1.249 

e* 

-.0184 

-.0071 

.0017 

-.1440 

-.1607 

-.1815 

-.1237 

-.1465 

-.2844 

.8000 

-.0189 

2.498 

-.047 

1.178 

1.131 


-.0219 

-.0114 

-.1584 

-.3280 


-.0042 

3.226 

-.014 

1.078 

1.064 

f 

-.0012 











-.1728 

-.1872 

-.1969 




lover surface 



f 

-.0210 

-.0130 

-.0031 

-.1329 

-.3298 








0.0125 

-0.0598 

8.263 

-0.494 

0.458 

-0,036 

6 

-*0186 

-.0164 

-.0035 

-.2093 

-.0818 

-.2911 

.0500 

-.0792 

4.603 

-.364 

-.314 

-.678 

s' 

-.0157 

-.0151 

-.0015 

-.2017 

-.2189 

.0060 

-.2129 

.1000 

-.0938 

3.304 

-.310 

-.541 

-.851 


-.0128 



-.2161 





-.1078 

2.413 

-.260 

-.676 

-.936 

h 

..0128 

-.0043 

-.2332 

.1268 

-.1064 

.4000 

-.1199 

1.914 

-.230 

-.754 

-.984 

h* 

-.0096 

-.0118 

-.0066 

-.2305 

-.2467 

.2747 

.0280 

.6000 

-.1073 

1.903 

-.203 

-.796 

-.999 






.8000 

-.0564 

2.335 

-.132 

-.834 

-.966 

1 

-.0071 

-.0092 

-.0122 

-.2449 

-.2642 

.4425 

.1783 

.9000 

-.0333 

3.080 

-.103 

-.845 

-.948 

J 

.0002 

-.0187 

-.0261 

-.2737 

-.2996 

.7982 

.4986 








.0115 

-.0362 

-.0323 

-.3025 

-.3233 

1 .U 87 

.7954 







D 

.0232 

-.0586 

-.0224 

-.3313 

-.3305 

1.3365 








■ 

.0281 

-.0687 

.0021 

-.3601 

-.3299 

1.4117 









a 

Velocities along the surface are considered posltlTe vhen directed from the 
trailing edge +c the leading edge on the lower surface, and from the 
loading edge to the trailing edge on the upper surface. , 
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TABLE VII 

CCMPUTATICK OF MAPPOTG FCNCTICIT CCHSTANTS 
FOE NACA kk 22 AIRFOIL LATTICE 


K 

®n 

®t 

S 

T 

U 




(deg) 

(deg) 




°^0 



- 7. 37 

181.62 

1.1307 

0.0157 

0.2879 

1 .4343 



-9-81 

182.38 

•9585 

.0152 

.2895 

1 .2632 


0.32 

-7.85 

181.83 

1.0920 

.0156 

.2886 

1.3962 


0.311 

-7.63 

181 .73 

1.1088 

.0157 

.2882 

1.4127 


0.308 

-7.57 

181.72 

1.1147 

.0157 

.2882 

1.4186 


0.309 

-7.59 

181.73 

1.1128 

•0157 

.2882 

1.4167 

s ^ 

0.3083 

-7.57 

181.72 

1.1144 

.0157 

.2882 

1.4183 h 


S « 


cos Xq 
2j(a 


logg 


(cosh K - COB 0^) (cosh K + cos 0j^) 
(cosh K + cos 0^) (cosh K - cos 0j^) 


(sin 0 J 1 - sin 0t)sinh. E 
slnh.^E + sin 6 ^^ sin 

TJ - 1* tan"l ("tag- % ~ tern. 0t)ts]ah K 

SjtcVQ tanh% + tan 9 ^ tan 0^ 


_ sin Xn ^ 

T ^ tan^ 


rtc 
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TABI^E TUI 


CO 

c» 


KEIATICH HEWEEH CIECIE M(XLE 6 (t-PLAITE) MD 
LOCATICW m HACA khl2 AIEFOIL HI C£SCA1® 


Upper surface 


Lower 

surface 

Point 


i/c 

e 

(^t^eg) 

Point 


x/o 

e 

( 4 eg) 

g 

-0.2911 

0.006 

1.4 

a 

0.5875 

1.000 

-178.3 

h 

-.1064 

.090 

17-0 

t 

.5400 

.920 

-167.3 

1 

.1783 

.270 

45’.7 

c 

.3106 

.707 

-135.0 

J 

.4966 

.501 

90.6 

a 

.0401 

.435 

-82.9 

k 

.7954 

.740 

145.2 

e 

-.2022 

.187 

- 35.1 

1 

1.0060 

.927 

168,2 

H 

-.3200 

.033 

-11.8 


TCT a rnTTr^r at* j u v ii ’ I' r t/vts v 

amxxLctiUf jwyxoux^i 

COMOrrEE POOR AKRCHAOnCS 
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fABLB IX 


caBSTAHTs car baca 4ina aihfoil latticb 



Source -TO(rtax aethod. 

Contour- 
Integral method 

1 

4'AW VUWVX V 

reference 8 

First 

ftppro 3 ciniEt-lo 3 i 

Seoeaod 

n nniv^-rt Tim.+.l mi 


0.006 

0.006 



dJ'^/cjQ 

-.538 

-.365 




-.324 

-.342 

- 0.346 



1.03 

.99 

.99 

1.00 

K 



.3083 

.3109 




ifii 7 ^ 

iftl 7Q 

^1 r 

dC2/acyj 



3.71 

3.71 

T} , deg 




-5.75 

-5.94 
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figure f . — flow S(r)guianfies m circle plane and corresponding velocity vectors 
m physical plane. 


i—»a 

dq 
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Fig. 2 
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Figure 2.- Velocity potential and stream function for a row 
of vortices of unit strength spaced at unit distance along 
the y-axis with the central vortex omitted. 
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Fig. 3 



Fig. 4 


NACA TN No. 1252 




— /o3a^- 




Figure. 4.- NACA 44 IF atrFot/ /n lattice 

arrangement. j3=^0° ^ C=/.03Z; 9.7 . 



— Isolated flow] Cf = /*67 
— /^f/jod of rsfierenoe 3Jc/-/.c?o 
Source-vortex, method'^ 

first approximation', q ~/.o3 
° Source-vortex method\ 

second, approx imatio/%q-a,93 
^ ContouMntegraf method\ I 
first approximation', q*o. 


tf. * tr / _ k t A jt ,n - / . _ ^ i JL _ ^ t 

ngure ~ yg/ocio/GS on /v/\or\ •^*tic. um on m mu/aifca rtatv onu 

lattice arrangement, <s*/.03B% oc*3.7^ , 




OQ 

o> 



•v-f 


44IZ « ^rsic^^rj: 

^ yari-tces alortg mean //>itf , 

chorf- readings taken. 


I \JUJ 

,/Z X 


li^~ 


/ A 

■ /-■ -:,■ 


. / > 

•■ ■ j — 


■ I" i- 


r .i ;■ 

« ■■ ;■■ 

n I 
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^. 4-0 







r t^ure ( 


/ — Soorce-uorier /n&ihod’^ 

yZ. ftrsf approjrimatior) 

y Source-vorie/ /neihod; 

— second approx/ma-h/on 

Contour-/ nfegrcil meihod 

i 

80 120 160 200 240 280 320 360 

Circie anqie , (p, deq NATiONAL advisory 

^ OOHMITTEE FW AEMMAUTICS 

A * 

The induced flow funcrion 'Xq agamsf arch ^ongie for N h Cf\ 
4412 Qirfo/I in lattice arrangement, /9 ~ 0; <r -z I.03Z ; 




Figure 6~ Typical curves for determination of ^and by contour- integral 
meth'od. These curves are for point g on NACA ^^12 airfoil m 
lattice arrangement. /3=o° <r= 1.032 ; cc^-5.7® 


Fig. 8 NACA TN No. 1252 


o 



-130 -140 -100 -60 -20 20 60 (OO 140 tSO 


Cl rdf angle j - 9 , deg 

Figure 9 Vdocti-g poienhot on unrt circle in plane ^ for NACA 4HZ 

airfoi] m iathce arrangement T 1032', cc^- 9'7®. 


Fig. 



r >gur& fO . Velocity distributions on NACA 

1.032 



(b) a^ = /a/®. 

44IB airfoil m lattice arrori^ement . 
; <^-- 0 °. 


* 


% B 


Fig. 10 N-AlCA TN No. 1252 
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Fig. 11 



O 

Figure ! / .—Denied oirfoif I ait ice . /9 - 45 ; 

cr--/. 006 ’^ 




along mean Une and locations at which chart readings were taren. 
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Fig. 12 NACA TN No. 1252 



Figure !3 Veloc/hes on cfer/yed pirFoif 
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O 40 60 }iO 160 100 140 160 610 660 


Circle angle, deg 

Figure 14. - The induced fkw fOnchon §(, against arde angle for -the denv^ed airfoil 
lattice, &‘-45\ (7=1006^ 
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iV? S~pfan9 for derived airfoil lattice 





